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Three-Body Model Analysis of Subbarrier a Transfer Reaction 

Tokuro FuKUi,^*) Kazuyuki Ogata, ^**) and Masanobu Yahiro^ 
^Department of Physics, Kyushu University, Fukuoka 812-8581, Japan 

Subbarrier a transfer reaction "C(^Li, d)^'^0(6.356 MeV, 1/2+) at 3.6 MeV is analyzed 
with a a + d+^^C three-body model, and the asymptotic normalization coefficient (ANC) 
for Qf+^'^C — > ^^0(6.356 MeV, 1/2+), which essentially determines the reaction rate of 
^''^C{a,n)^^0, is extracted. Breakup effects of ^Li in the initial channel and those of ^^O in 
the final channel are investigated with the continuum-discretized coupled-channels method 
(CDCC). The former is found to have a large back-coupling to the elastic channel, while 
the latter turns out significantly small. The transfer cross section calculated with Born 
' approximation to the transition operator, including breakup states of ^Li, gives {C^^q)^ = 

04 . 1.03 ± 0.29 fm^^. This result is consistent with the value obtained by the previous DWBA 

' calculation. 

. §1. Introduction 

Transfer reactions below Coulomb barrier energies are known to be a powerful 
technique to determine asymptotic properties of the overlap between the initial and 
final state wave functions, essentially free from uncertainties associated with optical 
potentials and structural complexity of wave functions in the nuclear interior regionP 
^ , Recently, subbarrier a transfer reactions have been used to indirectly measure cross 

Ch \ sections of a-induced reactions of astrophysical interest .'^''^ In Ref. [2]), Johnson 

and collaborators determined the reaction rate of ^^C(a,n)^^0 by measuring the 
i3c(6Li,d)i^O(6.356 MeV, 1/2+) reaction; for simplicity, we henceforth denote the 
final state of ^^O as ^^O*. The ^^C(a,n)^^0 reaction is considered to be important 
• as a neutron source for the slow neutron capture process (s-process) taken place in 

CN . the asymptotic giant branch (AGB) stars.l^ 

In the cross section formula, Eq. (1) of Ref.[2]), of the ^^C(a, n)^^0 reaction based 
■ on i?-matrix approach,'^ the asymptotic normalization coefficient (ANC) for a+^^C 

— > ^^O*, C^ipQ, is the only missing quantity. Throughout this study we consider 
the ANC with Coulomb-modification ,'2' i.e., a value divided by the Gamma function 
r{2 + rj), where rj is the Sommerfeld parameter for the a-^^C system. In Ref. [2]), 
k>( ■ the a transfer reaction ^^C(^Li, (i)^''0* was analyzed with DWBA, disregarding the 

5-1 ! breakup effects of ^Li and ^^O*, and {C^^i?^)^ = 0.89±0.23 fm^^ was extracted. The 

- ground state energy of ®Li is, however, just 1.47 MeV below the a + d threshold. Fur- 

thermore, the binding energy of ^^O*, i.e., ^^0(6.356 MeV, 1/2+), from the a-l-^^C 
threshold is only 3 keV. Therefore, to extract a reliable value of C^iPq, one should 
investigate how important ®Li and ^^O* breakup are in the a transfer reaction. 

The purpose of the present Letter is to analyze the ^^C(^Li, d)^^0* reaction at 
3.6 MeV (for the incident energy of ^Li) with the three-body (a + d+^^C) model and 
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to determine C'^iaJ^ accurately. Roles of ^Li breakup in the initial channel and ^''O 
breakup in the final channel are investigated with the continuum-discretized coupled- 
channels method (CDCC).'^''^ As shown in ^3.21 the former is found important as a 
large back-coupling to the elastic channel, while the latter is confirmed much less im- 
portant. CDCC was proposed and developed by Kyushu group and has been highly 
successful in quantitatively reproducing observables of reaction processes in which 
virtual or real breakup effects of the projectile are sig nificant.'3''3 CDCC treats con- 
tinuum states of the projectile nonperturbatively, with reasonable truncation and 
discretization, and thus can describe the breakup effects with very high accuracy. 
Note that theoretical foundation of CDCC was established in Refs. I10 p - tl2p . The 
transition from the ^Li+^'^C channel to the channel is described with Born 

approximation; the breakup states of ^Li are explicitly taken into account in the 
calculation of the transfer process. The ANC thus extracted is compared with the 
result of the previous DWBA analysis. 

This paper is constructed as follows. In ^ we formulate the three-body wave 
functions in the initial and final channels and the transfer cross section of the 
^^C{^Li,dy^O* reaction. Numerical setting is described in ^3.1[ Breakup effects 
of ^Li and ^^O are investigated in ^3.21 and the transfer cross section is analyzed 
and the ANC is extracted in ^3.31 In ^3.41 we see the convergence of the modelspace 
of CDCC, and in ^3.51 we discuss the present result in comparison with the previous 
DWBA result. Finally, we give a summary in ^ 

§2. Formulation 




Initial channel ~ " Final channel 



Fig. 1. Illustration of the three-body system in the initial and final channels. 

In the present calculation, we work with the a + d+^^C model shown in Fig. [TJ 
The transition matrix (T matrix) for the transfer reaction ^^C{^Li,d)^'^0* is given 
by 

Tj, = S'4{^^f^\VtA^^^^), (2-1) 

where ^^f"*"^ and iZ'j^ ^ are the three-body wave functions of the system in the initial 
and final channels, respectively, and Vtr is the transition operator of the transfer 

1/2 

process. We put a normalization constant Scxp in Tfi, physics meaning of which is 
discussed below. 

The three-body wave function '^^^^^ in the initial state satisfies the Schrodinger 
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equation 



0, 



(2-2) 



where E is the total energy of the system in the center-of-mass (cm.) frame and r 
(i?) is the coordinate of a (^Li) relative to d (^^C). The Hamiltonian Hi is given by 



rCoult 



(2-3) 



where Tr is the kinetic energy operator associated with R and hi is the internal 
Hamiltonian of ®Li. We use for the nuclear interaction between X and Y; each 
of X and Y represents a particle, i.e., d, a, or C (^'^C). Similarly, .RxY denotes the 
relative coordinate between X and Y. y^oui -g ^j^^ Coulomb interaction between ^Li 
and ^^C. Note that we neglect the Coulomb breakup of ^Li, which can be justified 
by the fact that the effective charge of the a + d system for electric dipole transition 
is almost zero. Furthermore, as shown in ^3.21 it is numerically confirmed that 
Coulomb breakup processes due to electric quadrupole and higher multipoles are 
negligibly small. 

As the partial wave ^i-jm 

, we adopt the following CDCC wave function: 



ymax -tmax 



CDCC 
i-JM 



j=o e=o L=\j~e\ 



R 



i%{r) ® i^YL{R) 



JAI 



(2-4) 



where J and M are the total angular momentum and its z-component, respectively, 
and £ (L) is the orbital angular momentum between a and d (®Li and ^^C). We 
disregard the intrinsic spin of each particle for simplicity. The radial part of the ^Li 
wave function is denoted by 4>j^e{r)/r, where j is the energy index; j = corresponds 
to the ground state and j / to discretized continuum states obtained by the 
momentum-bin discretization.El' The internal wave function ^j/^m given by 



(2-5) 



(2-6) 



satisfies 



Inserting Eqs. ()2-3|) and ()2-4p into Eq. (|2-2p and making use of Eq. ()2-6p . one 
obtains the following CDCC equation: 



X, 



i{R) = -Y,FcAR)xi{R), (2-7) 



2/x dR^ 2fi R^ 
where /i is the reduced mass of the ^Li-^^C system, Ejj = E — ej/, and 



Fcc'iR) 



JM 



y(N) ^(N) 



i' Yfi ® 



JM , 



r,R 

(2-8) 
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For simple notation, we denote the channel indices {j, i, L} as c. The CDCC equation 
is solved numerically up to = -Rmax and Xc is connected with the usual boundary 
condition 

<,J(m ^ I f^U,(^J/^)'^-o - VK^JjK-,Siy^l^{K,^,R) for E,,, > 
^ \ -Si^W_r,^^^L+i/2{-2iKj,,R) forE,-,<0 ' 

(2-9) 

where Kj^i = ^2[iEj^iijh^ ^Lr^ e (^L~»? f) incoming (outgoing) Coulomb wave 

function with the Sommerfeld parameter rjj^i, and W_^. ^^i_^_i/2 is the Whittaker 
function. The subscript of ^ and c represents the incident channel. With the S- 
matrix elements S^^^ in Eq. (j2-9p . one may obtain any physics quantities with the 
standard procedure except that one needs to make the discrete results smooth when 
breakup observables are investigated. 

Since the CDCC wave function tp'pDCC ^^^y^ regarded as, with very high accu- 
racy, an exact solution to Eq. ()2-2p in evaluation of T-matrix elements that contain 
a short range interaction, one may define Vtr by 

Vtr = Vad{r) + VdciRdc) + VaciRac) " Kux (2-10) 

with any choice of the auxiliary potential 14ux- In Eq. ()2-10p . Vad, Vdc, and VaC 
contain both nuclear and Coulomb parts. Note that laux determines the final state 
wave function iP'j. . In the present study, we adopt 

Kux = VdciRdc) + VaciRac) + (2-11) 

which trivially gives 

Kr = K!?(r). (2-12) 

The superscript (C) of V^d in Eq. (j2Tip represents the Coulomb part of the inter- 
action. We then have 

{Hf - E)^^j'-\Rac,Rdo) = (2-13) 

with 

Hf = Tr^o + VdciRdc) + KS?(r) + hf, (2-14) 

where Tr^q is the kinetic energy regarding -R^o and hf is the internal Hamiltonian 

of ^^O. Note that we here consider a Schrodinger equation for which is the 

time-reversal of l^'j \ 

One can easily obtain the form of based on CDCC, , just in 

the same way as in the initial channel, except that i) we should include Coulomb 
breakup of ^^O, ii) we have no nuclear part of Vad, and iii) the bound state of ^^O at 
6.356 MeV is a p-wave that generates both monopole and quadrupole interactions 
between d and ^^O; the latter causes also change in the angular momentum 

that is called reorientation. Note that Vdc in Eq. (j2-14p contains both nuclear and 
Coulomb parts, as mentioned above. 
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It is shown in ^3.21 that ^'^O breakup channels have very small (~ 5%) effects 
on the d-^'^0 elastic scattering. Furthermore, the quadrupole interaction is found 
negligibly small (see Fig. [2]) . Then we can approximate 

!^r''^-^-^o(r)Ci-^(i?.o)-^r^-\ (2-15) 

where ipo{r) is the relative wave function between a and ^^C in ^'^O*, and ^ (Rdo) is 
the distorted wave function obtained by the single-channel calculation, in which both 
the breakup channels and the aforementioned quadrupole interaction are switched 
off. 

In the calculation of Tjj, we make zero-range approximation; the strength Dj^i 
of the zero-range a-d interaction is given by 

^i,^ = / 4>li{r)V^Ti^)hA^)d^- (2-16) 

The finite-range correction to the zero-range calculation of Tfi is made with the 
standard prescription.l^ One may examine the validity of this approximation by the 
magnitude of the correction. We use calculated with CDCC, while wj""^^"^ of 

Eq. (|2-15p is adopted as \ in the evaluation of Tji . 

§3. Results and discussion 

3.1. Numerical input 

The a-d wave function in t^pDCC jg constructed by following Ref. [T3]) , except that 
we do not use the orthogonal condition model (OCM) but exclude Pauli's forbidden 
states by hand. We include i = 0, 1, and 2 states. As for the nuclear part of the a-d 
interaction for i = 0, we use 

K!S=o(^) = -105.5 exp[-(r/2.191)^] 46.22 exp[-(r/1.607)2]. (3-1) 

For e = 2, 

^i3=2(^) = -85.00 exp[-(r/2.377)2] + 30.00exp[-(r/1.852)2] (3-2) 

is adopted. We neglect the intrinsic spin 5 of d, and we have only one resonance state 
at 3.474 MeV (measured from the ground state energy) with a width of 0.45 MeV. It 
is found that, however, if we include S and a spin-orbit interaction that reproduces 
the l"*", , and 3"'' resonance states, the resulting value of the ANC shown below 
changes by only 0.2%. Thus, the separation of the i = 2 resonance state to 1"*", 2+, 
and 3"^ resonance states by the spin-orbit interaction plays no role in the present 
subbarrier a transfer reaction. For £ = 1, we adopPf 

= -74.19exp[-(r/2.236)2], (3-3) 

which is used also for i > 2 when we check the convergence of CDCC calculation with 
respect to ^max (see ^3.4|) . The Coulomb interaction between a and d is evaluated 
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by assuming a uniformly charged sphere with the charge radius Rc of 3.0 fm; see 
Eq. ([33]) below. 

We take the maximum value /cmax (^max) of the relative wave number k (coordi- 
nate r) between a and d to be 2.0 fm~^ (60 fm); the maximum relative energy Cmax 
is 62.4 MeV. We use jmax = 100 for each of the £ = 0, 1, and 2 states and the width 
Ak of the momentum bin is thus 0.02 fm~^. The number of channels, Nch, in the 
CDCC equation ()2-7p is 601. When we see the effects of Coulomb breakup in Fig. [21 
we take r^^x = 300 fm. 

As for the interactions of the a-^^C and d-^^C systems, we use the parameters 
shown in Table HI The standard Woods-Saxon form is adopted: 

V{x) = -VoMx) - iWoMx) + Vc{x), (3-4) 

where fy{x) = (1 + exp[(3; — Ry) / ay])~^ and fw{x) = (1 + exp[(x — i?w)/aw])~^- 
The Coulomb interaction Vc{x) is given by 



Vc{x) 



3 - -2- ) x<Rc 



2Rc V Rc 



(3-5) 

x > Rc 



where Z1Z2 is the product of the atomic numbers of the interacting particles. These 
parameters are used in the calculation of both initial and final state wave functions. 
The parameter set for the d-^^C system is determined to reproduce the elastic scat- 
tering cross section obtained with the parameters in Ref.[2]) that contains a spin-orbit 
part. We determine Vq for the a-^'^C system to reproduce Eq assuming that the or- 
bital angular momentum is 1 and the number of forbidden states is 2. Note that 
we use Eq. ()3-5p with Rc = 2.94 fm for the ^Li-^^C Coulomb interaction unless we 
include Coulomb breakup of ^Li. 



Table I. Potential parameters used in the present calculation. 



System 


Vo 


Rv 


ay 


Wo 


Rw 


aw 


Rc 




(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) 


(fm) 


(fm) 




69.30 


2.939 


0.670 








2.969 




73.05 


3.128 


0.780 


10.50 


2.986 


0.800 


2.969 



In the calculation of ^p^^^, we use i?max = 15 fm and Jmax = 7. Note that 
we explicitly include closed channels, in which Ej^i < 0, in CDCC calculations. In 
the evaluation of Tjj, we set the maximum value of Rdc to be 30 fm; we use the 
asymptotic form of Xc ) Eq. (|2-9p . to obtain ^pj^^*^ for R > 15 fm. When we include 
Coulomb breakup, we set -Rmax to 200 fm. 

For the final channel, the relative energy between a and ^'^C in the 1/2"'" state at 
6.356 MeV is Eq = —3 keV from the a-^^C threshold. In the calculation of 1^^^^^^ ^ ^ 
we include the p-wave bound state and the s-, p-, d-continua of the a+^^C system 
up to the relative momentum of 1.2 fm~^ (relative energy of 39.6 MeV) with the 
momentum bin with a common width 0.06 fm^^. The maximum values of RaC and 
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Rdo are both set to 100 fm, and we put Jmax = 10. We include all closed channels 
in the CDCC calculations as in the initial channel. 

3.2. Breakup effects of^Li and O 




30 60 90 120 150 180 30 60 90 120 150 180 



e (deg) e (deg) 

Fig. 2. (color online) Elastic cross sections of ^Li-^''C at 3.6 MeV (left panel) and at 
1.1 MeV (right panel). In each panel, the dashed line shows the result of CDCC and the solid 
line is the result without breakup channels. The dotted line in the left panel is the result of 
CDCC with both nuclear and Coulomb breakup, while that in the right panel shows the result 
without breakup including only the monopole interaction between d and ^'^O. 

Figure [2] shows the elastic cross sections of ^Li-^'^C at 3.6 MeV (left panel) and 
^_17q* j MeV (right panel) corresponding to the initial and final channels, 
respectively, of the ^^C(^Li, (i)^'^O* reaction. In each panel, the dashed line shows 
the result of CDCC and the solid line is the result without breakup channels. One 
sees from the left panel significant breakup effects on the elastic cross section, i.e., 
a large back-coupling to the elastic channel. Another finding is the inclusion of 
Coulomb breakup (the dotted line in the left panel) little affects the cross section. 
One can thus infer that nuclear breakup plays important roles in the ^^C(^Li, d)^^0* 
reaction, and conclude that neglect of the Coulomb breakup in the calculation of 
^PFm'^ is justified. On the other hand, in the final channel, effects of nuclear and 
Coulomb breakup are found very small as shown in the right panel; they change the 
cross section for 6 > 60° by 5% at most. We further investigate the breakup effects 
on the d-^^O* wave function in the elastic channel. The absolute value (argument) 
of the wave function for J = at Rdo = 10 fm, which is found to have the main 
contribution to the transfer amplitude, is 0.982 and 0.956 (278.8° and 276.7°) when 
the breakup states of ^^O are included and neglected, respectively; the breakup 
effects are about 3%. Therefore, we can disregard the breakup channels of the d-^^O 
system in the calculation of Tjj with the error of 5% at most. The dotted line in 
the right panel shows the result with neglecting both the breakup channels and the 
quadrupole interaction between d and ^^O, which is almost identical to the solid 
line. Thus, one can use Eq. (|2-15p in the calculation of the final state wave function; 
we estimate the error due to this approximation to be 5% as mentioned above. It 
should be noted that breakup cross sections in the initial and final channels are both 
found smaller than the nuclear part of the elastic cross section by about four orders 
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of magnitude. 

The very small breakup effects in the final channel are because the incoming 
energy of d is suitably below the Coulomb barrier, and the interaction that causes 
breakup in Eq. ()2T4p is significantly weaker than in Eq. (j2-3p : note that y^\r) is 
defined as Vtr and does not appear in Eq. ()2T4p . 

3.3. Transfer cross section and ANC 





e (deg) 

Fig. 3. (color online) Cross section of the transfer reaction ^^C(^Li, d)^^0* at 3.6 MeV. The solid 
line is the result of calculation with Soxp = 1- The dashed line is the result of the fit to the 
experimental data taken from Ref. [2|. 



We show in Fig. [3] the cross section of the transfer reaction ^^C(^Li, (i)^''0* at 
3.6 MeV as a function of the outgoing angle of d in the cm. frame. The solid line 



represents the result with 5, 



exp 



1 and the dashed line shows the result of the x 



fit to the experimental data.'^J The resulting value of S^m is 0.357. Note that 5, 



-•cxp 

cannot be regarded as a spectroscopic factor. Indeed, Scxp has strong dependence 
on the model wave function of the a-^^C system; typically, it varies by a factor of 2 
with changing the geometric parameters of V^^^ by 30%. This clearly shows that it 
is not feasible to determine S'exp from the present analysis of the experimental data. 
On the other hand, the ANC C^^^q given by 



n^'^O* _ qI/2 ^sp 



(3-6) 



with the single particle ANC C^^^o,^ of the a-^^C wave function, is robust against 
changes in the potential parameters. This shows that the reaction process considered 
is peripheral with respect to Rac^ only the tail of the a-^^C wave function 
contributes to the transition amplitude. Note that C'^isq is defined by 



/^sp 



W^-f),3/2(2Ato^«c)r(2 + j7) 



at Rac ^ -Rn, 



(3-7) 



where <^o is the radial part of ipQ, rj is the Sommerfeld parameter of the a-^^C system. 
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Ko = ^— 2/iQi3c £o/h with /iQ,i3c the reduced mass of a and ^"^C, F is the Gamma 
function, and i^N represents the range of . 

The value of (^^1?^)^ extracted by the present calculation is 1.03 fm^^. We then 
evaluate the uncertainty of (C^^^q)^ associated with the a-^^'^C and d-^^C potential 
parameters shown in Table |T]by changing each value by 30%. Note that Vq for a-^^C 
has a constraint that it must reproduce Eq. The uncertainty is found to be 22%. 
We take into account also the uncertainty due to the use of Eq. ()2-15p (5%) and 
that coming from the zero-range approximation to (8%), and conclude that 

the theoretical uncertainty is totally 24%. Including the ambiguity of experimental 
informatiorP together, we finally obtain 



1.03 ± 0.25 (theor) ± 0.15 (expt). 



(3-8) 



where (theor) and (expt) respectively represent theoretical and experimental uncer- 
tainties. 

3.4. Convergence of the CDCC wave function in the initial channel 



G 



0.1 



0.05 



(0) 

0.26 (1.05) 

■ 0.40 (2.50) 

■ 1.6 (39.9) 

■ 2.0 (62.4) 

■ 2.4 (89.9) 




0.1 



0.05 



s,p, d; Ayt=0.02 
s, p, d; Ayt^O.Ol 
s,p,d, f, g,h; tJi' 
without breakup 
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Fig. 4. (color online) The dependence the cross section for ^■^C(''Li, d)^^0* at 3.6 MeV on the 
modelspace of CDCC for <Z'P°'^*^. In the left panel, fcmax is varied. The values of fcmax are 
shown in unit of fm~^ and the corresponding values of tmax are given in the parentheses in unit 
of MeV. Here, the = 0, 1, and 2 breakup continua are taken with Ak = 0.02 fm~^. In the 
right panel, the dashed line stands for the result of the ^ = 0, 1, and 2 breakup continua with 
fcmax = 2.0 fm~^ and Ak = 0.01 fm~^. The dotted line shows the result of the < £ < 5 
continua with fcmax = 2.0 fm~^ and Ak — 0.02 fm~^. The thick solid line (the solid line) in the 
left (right) panel is the result of the ^ = 0, 1, and 2 breakup continua with fcmax = 2.0 fm^^ and 
Ak = 0.02 fm"^ and the same as the solid line in Fig. [31 The result without breakup channels 
is also shown by the dash-dotted line in each panel. 



FigureSlshows the dependence of the cross section for ^^C(^Li, d)^'^0* at 3.6 MeV 
on the modelspace of CDCC for ip'pDCC^ ^^iq left panel, we show the convergence 
of the cross section with respect to increasing /cmaxj where the i = 0, 1, and 2 breakup 
continua are taken with Ak = 0.02 fm~^. One can see that the convergence is very 
slow and obtained at /cmax = 2.0 fm~^. In usual CDCC calculation, one takes only 
the open channels, i.e., channels with Ej^i > 0. The result thus obtained (the thin 



10 



T. Fukui, K. Ogata, and M. Yahiro 



solid line) is, however, sizably different from the converged one (the thick dotted 
line), at backward angles in particular. Thus, inclusion of the breakup channels is 
important. 

In the right panel of Fig. [U the dashed line is the result including the i = 0, 1, 
and 2 breakup continua with zA^ = 0.01 fm ^ and /cmax = 

2.0 fm-i (iVch = 1201), 

and the dotted line is the result including the ^ = 0, 1, 2, 3, 4, and 5 breakup 
continua with = 0.02 fm ^ and A;]-[iax = 

2.0 fm-i (A^ch = 2101). The dashed and 
solid lines both agree well with the solid line, which is the same as in Fig. [3l In fact, 
the resulting values of (C„i?c)^ differ from each other by less than 1%. Thus, the 
modelspace used in the solid line of Fig. [3] gives good convergence of the calculated 
cross section, hence C i9^. 

3.5. Discussion on the comparison with the previous DWBA analysis 




I I 

30 60 90 120 150 180 

e (deg) 

Fig. 5. (color online) ®Li-breakup effect on the transfer reaction ^^C(®Li, ci)^^0* at 3.6 MeV. 
The solid and dash-dotted lines show the results of calculations with and without ''Li breakup 
channels, respectively. The dashed line stands for the result of the elastic transfer process only. 
The solid and dash-dotted lines are respectively the same as those in the right panel of Fig. [l] 

Figure [5] shows the ^Li-breakup effect on the transfer reaction. The solid and 
dash-dotted lines are the results of the calculation with and without the breakup 
channels, respectively; they are shown also in the right panel of Fig. HI The two 
results largely deviate from each other, indicating that the breakup effect is impor- 
tant, as mentioned in ^3.4[ This does not necessarily mean, however, inadequacy of 
DWBA, as discussed below. The transition matrix elements of the transfer reaction 
can be separated into two parts, 




with 

rjf = (^}-)|T4,|a^W), (3-10) 
T)^=-) = (u.(-)|y,,|^W), (3-11) 
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where '^^^ and iP'j^j, are the elastic and breakup parts of the CDCC wave function 
^CDCC^ The transition matrix T^f^ describes the transfer reaction from the elastic 
channel, i.e., the elastic transfer process, which includes the back-coupling effect of 
the breakup channels to the elastic channel. On the other hand, T^Y^ describes 
the transfer reaction from ^Li breakup channels, i.e., the breakup transfer process. 
Thus, there are two kinds of breakup effects on the transfer reaction; one is the back- 
coupling effect in the elastic transfer process and the other is the presence of the 
breakup transfer process. The dashed line is a result of the elastic transfer transition 
only. The result agrees with the solid line, indicating that the breakup transfer 
transition is much smaller than the elastic transfer one. This is consistent with the 
small breakup cross section of ^Li by ^^C as mentioned in ^3.21 Hence, only the back- 
coupling effect is important in the present subbarrier transfer reaction. In DWBA, 
the back-coupling effect is expected to be included by using the ^Li optical potential, 
which describes the elastic scattering by definition, as the distorting potential. The 
ANC, (C^i?c)2, extracted in the preceding DWBA calculatiorP is 0.89 ±0.23 fm'^ 
This value agrees well with the present result Eq. (j3-8p within the uncertainties. 



§4. Summary 

In summary, we analyze the ^^C(^Li, d)^ ''0(6.356 MeV, 1/2'^) reaction at 3.6 MeV 
by the three-body (a + d+^^C) model. The breakup effects of ^Li and ^^O are in- 
vestigated by CDCC. Those of ®Li are found important as a large back-coupling to 
the elastic channel, while those of ^^O turns out negligible with an error of 5%. The 
transfer cross section is calculated with Born approximation to the transition inter- 
action, and including only the breakup of ^Li. The ANC extracted by the three-body 
reaction model is (C^i^c)^ = 1-03 ± 0.25 (theor) zt 0.15 (expt). The back-coupling 
effect of ^Li breakup on the transfer reaction is large, while the breakup transfer 
transition is negligible compared with the elastic transfer transition. The preceding 
DWBA calculation implicitly treated the back-coupling effect by using a ^Li optical 
potential that described the elastic scattering as the distorting potential. The value of 
{C^i^c)"^ extracted by DWBA is 0.89ib0.23 fm~^, which is consistent with the present 
value within the uncertainties. It can be conjectured that in the DWBA calculation, 
the aforementioned back-coupling effect in the initial channel was properly included. 
However, this will not always be the case, since the optical potential is determined 
phenomenologically. Furthermore, breakup transfer processes may be important in 
other subbarrier a transfer reactions. The present three-body approach, therefore, 
should be applied systematically to these reactions. From theoretical point of view, 
inclusion of CDCC wave functions in both initial and final channels will be an im- 
portant subject; to achieve this, one should treat, in principle, very large coordinate 
space in the calculation of the T matrix, since there is no damping in the overlap 
kernel. It will be interesting to use four-body CDCd^ based on a p + n + a+^'^C 
model to obtain the wave function in the initial channel. At this stage, however, the 
modelspace required is too large for four-body CDCC to be applied. 
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